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This paper describes and discusses a general method of solutions for treating effects of 
Doppler-broadened Rayleigh backscattering in H 2 0 DIAL measurements. Errors in vertical 
DIAL measurements caused by this laser line broadening effect can be very large 
(Ansmann, 1) and, therefore, this effect has to be accounted for accurately. 

To analyze and correct effects of Doppler-broadened Rayleigh backscattering in DIAL 
experiments, a generalized DIAL approximation was derived starting from a lidar 
equation, which includes Doppler broadening: 
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where N(R) is the water vapor number density averaged over the range cell DR centered at 
R = (Ri + R 2 )/2 and Pi(R) and P 2 (R) are the on- and off-line laser signals backscattered from 
the range cell at R, respectively. The absorption cross section for wave number v, (i = 1 : on- 
line, i = 2: off-line) averaged over the range cell DR centered at R is given by 
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and the Doppler broadening correction term D, for wave number v, is given by 
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where 6 (Rj) is the total backscattering coefficient for range R Jf Bm(Rj) and 8 r(Rj) are the 
Mie and Rayleigh backscattering coefficients for range R J# h,(v) is the normalized laser line 
shape for wave number v u and g,(v, R^ is the convolution of the incident laser profile h ( (v) 
and the intensity distribution of an incident monochromatic frequency backscattered by 
air molecules in the range cell at Rj. 

The calculation of water vapor density N(R) with Eqs. ( 1 )-(3) includes the correction of 
Doppler broadened Rayleigh backscattering effects. For using this calculation scheme, laser 
and H 2 0 absorption line parameters, temperature and pressure have to be known as usual 
for H 2 O DIAL retrieval. In addition, for the correction of effects of Doppler broadening, 
backscattering properties of molecules and aerosols have to be known. The backscattering 
properties of molecules can be determined from measured temperature and pressure pro- 



files or approximated from model atmosphere data. The aerosol backscattering properties 
can be obtained from range normalized off-line signals P 2 (R)-R 2 3 4 5 using a numerical 
integration scheme (Fernald, 2; Sasano et at., 3). Here, in addition, Rayleigh backscattering 
coefficients, aerosol extinction/backscattering ratio and aerosol backscattering coefficient 
at calibration range R c have to be known. To yield generally stable solutions of this 
integration scheme for 8 m(R) with respect to uncertainties in the signals and the other in- 
put parameters, backward integration mode with boundary value Gm(Rc) f° r maximum 
range R c = R m ax should be used (Klett, 4). For the evaluation of the integral in Eq. (3) giving 
the water vapor optical thickness up to the range Ri, values for N(r) calculated in the 
preceding steps Ro<R<Ri can be used. The off-line Doppler broadening correction term 
D 2 can be neglected, if off-line absorption is approximately constant over the transmitted 
spectrum. 

To evaluate the accuracy of H 2 O DIAL measurements using Eqs. (1)-(3), computer 
simulations were performed using Gaussian laser line shape with HWHM of 0.0125 cm 1 , 
Voigt-H20 absorption lines with Lorentz halfwidth of 0.09 cm -1 , profiles of temperature, 
pressure and H 2 O density (5.9 g m-3 at ground level) for standard atmospheric conditions, 
aerosol scattering coefficients for clear atmosphere conditions (McClatchey, 5, ground 
level visibility of 23 km) with additional layers of enhanced aerosol backscattering, and 
realistic lidar and receiver parameters. Lidar signals needed for the determination of H 2 0 
and aerosol backscattering were calculated with a lidar equation considering Doppler 
broadened Rayleigh backscattering. H 2 O-DIAL number densities were calculated with Eqs. 
(1)-(3). To analyze only Doppler broadening effects, no errors in the profile ln(P 1 (R)/P 2 (R)) 
used for solving the DIAL approximation (1) were assumed. The difference between the 
retrieved and the true or input H 2 0 is shown in the figures. 

In Fig. 1 examples of error profiles to be expected in ground based H 2 0 DIAL 
experiments due to realistic uncertainties in input parameters for backscattering retrieval 
are shown. It can be seen that correction of Doppler broadening effects is generally 
necessary. The accuracy of this correction mainly depends on the accurate determination 
of aerosol backscattering coefficients. Analysis of the influence of each input parameter 
for aerosol backscattering retrieval shows that it is especially important to have a proper 
estimate of the boundary value (3 m(Rc) R c = R max as a starting value for iterative aerosol 
backscattering retrieval. Error peaks at the upper layer boundaries in Fig. 1 are mainly due 
to errors in 6 M (R C ) The influence of errors in backscattering calculation is most important, 
when Mie and Rayleigh backscattering coefficients are of the same order of magnitude 
and, in addition, layers with steep gradients of aerosol backscattering are present. This is 
the case in the figures for the layer regions in the middle and upper troposphere. The 
relative water vapor error due to uncertainties in Doppler broadening correction is 
proportional to DR- 1 so that the shown effects decrease with increasing DR. 

In conclusion, correction of Doppler broadened Rayleigh backscattering is possible 
with good accuracy in most cases of tropospheric H 2 0 DIAL measurements, but great care 
has to be taken when layers with steep gradients of Mie backscattering like clouds or 
inversion layers are present. 
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E'Q 1 : Errors in H 2 O DIAL measurements due to Doppler-broadened Rayleigh 

backscattering. Dashed lines: without correction of Doppler broadening; solid lines: 
correction applied. Errors in input parameters assumed for retrieval: 3K in temperature, -10 
hPa in pressure; aerosol extinction/backscattering ratio assumed for retrieval: 60 sr, true 
value: 80 sr; Mie/Rayieigh backscattering ratio at Rma* assumed for retrieval: 0, true value: 
0.18. Aerosol layers with an aerosol concentration a factor 2 (o) and 5 (•) higher than the 
aerosol concentration of the clear air model between 0.5-1 km (DR = 50 m), 1-2 km 
(DR = 200 m), 2-4 km (DR = 500 m), and 4-8 km (DR = 1000 m). 
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